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ON  THE  BAROTROPIC  MODEL  OF  THE  OCEAN  CIRCULATION 


V.  BARCILONt,  P.  CONSTANTINtAND  E.S.  TITIt 
Abstract.  We  examine  the  mathematical  properties  of  the  solutions  of  a  barotropic, 
wind  driven  ocean  with  bottom  friction  on  both  a  0—  and  and  /—plane  .  Except  for  small 
Rossby  numbers,  the  uniqueness  of  the  solutions  of  the  corresponding  partial  differential 
equations  dependent  on  an  a  priori  bound  for  the  gradient  of  the  vorticity.  For  the  /—plane, 
two  drivings  are  considered  which  give  rise  to  explicit,  global  unique  solutions.  For  large 
Rossby  numbers,  a  novel  nonlocal,  nonlinear  boundary  value  problem,  which  does  depend 
on  the  ^-effect,  is  obtained  for  the  circulation.  , 

1.  Introduction.  This  paper  is  concerned  with  the  question  of  whether  ocean  cir¬ 
culation  models  have  unique  steady  solutions.  This  question  is  very  important  for  under¬ 
standing  the  long-term  behaviour  of  these  models.  More  specifically,  we  shall  consider  this 
question  for  the  simplest  such  model,  namely  that  of  a  homogeneous  wind-driven  ocean, 
with  bottom  friction  and  no  topography. 

For  the  readers  unfamiliar  with  the  oceanographic  jargon,  the  above  ocean  circulation 
model  assumes  a  simplified  dynamics  embodied  in  the  following  set  of  equations: 

I  ^  j  I  -I-  uuj  -  fv  =  -fypx  -  Ku-  Ti,  'j 
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These  equations  govern  the  two-dimensional  motion  of  a  homogeneous  fluid  viewed  from 


a  rotating  frame  of  reference.  The  usual  oceanographic  notations  are  used,  namely  u 


and  V  represent  the  velocity  components  in  the  x  and  y  directions  which  point  eastward 


and  northward  respectively;  p  is  the  density;  p  the  dynamical  pressure,  and  (Ti,r2)  the 


wind  stresses  in  the  eastern  and  northern  direction.  The  Coriolis  parameter  /  is  equal 


to  /o  +  f3y.  Finally,  K  is  a.  Rayleigh  friction  coefficient.  Since  there  are  many  instances 


where  a  detailed  analysis  of  the  effects  of  the  bottom  viscous  boundary  later  give  rise 


to  Rayleigh-lilce  friction  terms,  this  type  of  dissipation  is  referred  to  as  “bottom  friction”. 


Thus,  these  equations  include  all  the  physical  processes  which  are  believed  to  be  important, 


namely  inertial  forces,  Coriolis  force,  dissipation  and  wind  forcing.  In  spite  of  the  fact  that 


these  equations  constitute  an  enormous  oversimplification  of  the  oceanographic  reality. 


mathematiciUly  they  are  stiU  very  difficalt  to  deal  with  because  of  their  nonlinearity. 


In  order  to  avail  ourselves  of  the  a-priori  estimates  for  the  strength  of  the  circulation, 


for  the  magnitude  of  the  currents,  and  for  the  vorticity,  which  we  previously  derived 


(Barcilon  et  eil.,1987),  we  switch  to  the  foUoAving  dimensionless  formulation  of  the  problem: 


VV  =  u; 


in  n 


(1.2a) 


=  0  on 
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The  notations  are  identical  to  those  of  Veronis  ( 1966a, b)  whose  early  numerical  study  of 


this  problem  has  been  very  inspirational  to  us.  The  Rossby  number  R  represents  a  relative  ; 


measure  of  the  inertial  effects;  e,  which  is  akin  to  the  Ekman  number,  represents  a  relative  y  CodoJ 

and/or 

2  ;  xiit  i  Special 
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measure  of  frictioncil  effects.  The  circulation  is  described  by  the  streamfunction  V’  and 
the  vorticity  u;.  The  only  differences  with  Veronis’s  work  are  (i)  that  the  curl  of  the  wind 
stress,  which  we  denote  by  F,  is  taken  to  be  a  general  function  of  x  cind  y  rather  than 
sin  X  sin  y  and  (ii)  the  beisin,  denoted  here  by  Cl  with  co2istline  dCl,  is  taken  to  have  an 
arbitrary  shape  rather  the  square  0  <  i,  y  <  ir  . 

2.  Uniqueness:  We  show  in  this  section  how  the  uniqueness  hinges  crucially  upon 
whether  or  not  the  gradient  of  the  vorticity  is  bounded,  or  in  mathematical  terms,  upon 
whether  we  can  get  an  a  priori  estimate  of  llVu;|loo-  As  usual,  |1  ||oo  stands  for  the 
norm  of  a  function,  namely 


Halloo  =  sup|/l| 

o 

Let  and  be  two  solutions  of  (1.2).  If  we  write: 

rl;'  —  il;  +  (i>  (2.1a) 


and 


u'  =  u;  +  T} 


(2.16) 


then  clearly 


T])  +  RJ{<i>,uj)  +  RJ((i>,  77)  +  +  e77  =  0 


vV  =  V 


(2.2) 


where  J{f,9)  stands  for  the  Jacobian  -  f^gx-  Multiplying  (2.2a)  by  and 


integrating  over  the  basin,  we  deduce  that: 


since 


R  I  (i>J{ip,T])dA-r  f  <f>(j>xdA  +  e  f  (l)V^<f>dA  =  0 

Jsi  Ja  Ja 

[  <pJ(4>,h)dA=  f  {{h^hy),  +  {-h^hMdA 
Jn  Ja  ^  ^ 

—  I  {nihy  —  Tiihx}  ds 

Jdil 


(2.3) 


=  0 

in  view  of  the  fact  that  <()  vanishes  on  the  boundary.  In  the  above  formulas,  h  is  an  arbitrary 
function,  n  =  (ni,n2)  is  the  outer  unit  normal  and  s  is  the  arc  length  along  the  coast. 
Making  use  of  the  divergence  theorem  on  the  last  two  terms,  (2.3)  can  be  written  as 


e||V0|l' 


-L 


•(77Q)d.4 


(2.4) 


where 

Q  =  (u,u)  =  (-V’y.V’i)  (2.5) 

is  the  velocity  field  associated  with  V*-  The  notation  ||  ||  is  the  customary  shorthand 

for  the  X2(fl)  norm  of  a  function,  namely 


||4f  =  !  \h\'‘dA 
Ja 

The  term  on  the  left  hand  side  of  (2.4)  can  be  transformed  as  follows.  First,  we  use  the 
divergence  theorem  once  again  to  get 


/  fl>V  (T7Q)d>l  =  -  f  rjq-V((>dA 
Jn  Ja 


then,  taking  absolute  values  of  both  sides  we  immediately  see  that 


\ri\\^<f>\dA 

In  deriving  this  formula  we  have  used  the  obvious  fact  that  the  modulus  of  Q  and  that  of 
Vt/;  are  identical.  The  last  step  consists  is  using  Cauchy’s  inequality  to  overestimate  the 
right  hand  side  in  the  l<ist  equation.  This  implies  that 

\  [  ^^■{nQ)dA\<\\vnM\  liv^ll 

Jn 

Substituting  this  estimate  in  (2.4),  we  conclude  that 


f  <?iV-(;?Q)d41<||V0||oo  f 
Ja,  Ja. 


c||V<i||<iJ||ViMUI|r/|| 


(2.6) 


Similarly,  multiplying  the  first  equation  in  (2.2)  by  rj  and  integrating,  we  see  that 


because 


R  f  TjJ{<p,u;)dA  +  f  T}(px  dA  +  €\\t]\\'^  —  0 
Ja  Ja 

J^TiJ{h,ri)dA  =  “  h,{^T]^)x}dA 

Jda  ^  ^ 

=  0 


if  k  stands  for  either  0  or  V>  which  both  vanish  on  the  coast.  Applying  as  before  Cauchy's 
inequality,  we  see  that. 


e|hl|<(i2||Vu;||„  +  l)||V^i| 


(2.7) 


Combining  the  two  inequalities,  we  get: 


ilV^ll  <  ^||VvlU(i  +  i2||Vu;l|„)iiV.^ll  (2.8) 

Thus,  if 

|iW|„(l  +  J!||VHU)<l  (2.9) 

V and  hence  0  and  r}  are  identically  zero,  i.e.  the  solution  is  unique.  We  should  recall  at 
this  stage  that  H-^Hoo  contains  e.  Unfortunaly,  we  do  not  have  an  expression  for  |lVu;||oo. 
Therefore,  (2.9)  need  not  coincide  with  the  regime  dominated  by  dissipation. 

The  uniqueness  in  this  viscously  dominated  regime,  can  be  deduced  by  writing  (1.2) 
as  an  integral  equation,  namely; 

0=  /*  FGdA-hR  f  (2.10) 

Ja  Ja 

where  G  is  the  Green  function 

G,  =6(r-r')  (2.11) 

For  small  R,  we  can  show  by  classical  arguments  that  a  power  series  in  R  converges. 

In  view  of  the  non  linearity  of  the  problem,  it  would  not  be  surprising  if  the  viscous, 
Stommel  like  solution  were  to  bifurcate  for  a  certain  value  of  R.  In  the  next  section,  we 
shall  consider  the  /—plane  problem,  which  though  very  different  oceanographically  from 
the  — plane  problem,  is  in  fact  rather  similar  mathematically.  In  particular,  progress  on 
the  question  of  uniqueness  is  also  hampered  by  the  lack  of  an  a  priori  bound  on  the 
gradient  of  vorticity.  Therefore,  it  is  surprising  to  see  meaningful  wind  stresses  yielding 
unique  circulations  for  all  vedues  of  R. 


3.  The  /—plane  problem:  The  analog  of  equations  (1.2)  are 


with 


R{i}xUjy  —  ifytjJx)  +  euj  =  F 


=  (jj 


in  U 


(3.1a) 


•0  =  0  on  dVl 


(3.16) 


As  mentioned  previously,  the  uniqueness  derivation  leads  to  an  inequality  similar  to  (2.8), 
namely 

>  i?'||VT6||cc  ||Vu;|U  (3.2) 


The  reason  for  considering  the  /—plane  problem  is  that  because  of  its  simplicity,  we  are 
able  to  exibit  two  solutions  valid  for  <dl  R's.  These  solutions  are  the  responses  to  the 
forcings  /  =  —  1,  and  F  =  —ui,  where  ui  is  the  first  eigenfunction  of  the  Laplacian  for 
the  basin.  Incidentally,  it  is  ironic  that  Veronis  used  the  forcing  F  =  sin  x  sin  y,  which  is 
in  fact  the  first  eigenfunction  for  the  square  basin  0  <  i,  y  <  ir !  However,  his  analysis  was 
for  the  /3— plane. 

We  start  with  the  case  F  =  —1,  i.e.  of  a  constant  wind  curl.  If  we  set 


a;  =  -e-'  (3.3) 

then  the  vorticity  equation  is  identically  satisfied.  The  streamfnnction  is  found  by  solving 

=  -c-^  (3.4) 


subject  to  the  boundary  condition  along  the  coast.  For  this  particular  solution  {|Va;|{oo  is 
identically  equal  to  zero.  Hence,  this  is  the  unique  solution  valid  for  all  Rossby  numbers. 


Incidentally,  the  lack  of  dependence  on  the  Rossby  number  should  be  noted  and  it  should 
be  recalled  that  cill  our  estimates  on  the  norms  of  ij)  and  uj  obtained  in  Barcilon  et  al. 
(1987)  were  independent  of  R. 

The  case  F  =  —  ui,  where 

V^Un  +  AnUn  =  0,  in  Q 


(3.5) 


■Un  =  0  on  dil 


is  more  subtle.  Once  again  we  can  check  that  the  following  expressions  satisfy  the  equations 
(3.1). 


i’  =  - 


Unix,y) 
eXr.  ' 


j}  = 


We  then  go  through  the  uniqueness  procedure  and  write 


,  ^ix,y)  ,  ^  ' 

i)  = - -j —  +  0, 


e.\„ 

«n(a;,y) 


+  n 


By  means  of  some  simple  manipulations  we  see  that: 


R 


and 


—  (Un,J(<^,7?))  =  0, 


--(«n.^('^,»7))  +  e|kll^  =  0 


This  implies  that 


(3.6) 


(3.7) 


Il0il  =  A, 


(3.8) 


i.e.  that  0  is  a  linear  combination  of  the  eigenfunctions  associated  with  A„.  Thus,  the 
multiplicity  of  solutions  associated  with  the  driving  F  =  tt,  is  equal  to  the  multiplicity  of 


the  eigenvalue  A^.  Interestingly,  the  first  eigenvalue  has  multipUcity  one  (Protter  1987). 
Thus,  for  F  =  -ui,  the  solution  exhibited  above  is  unique  for  all  R's. 


4.  Large  Rossby  number  asymptotics;  If  we  multiply  the  governing  equation 

(1.2)  by  (i/’xWj,  -tpyf^x)  which  we  denote  for  brevity  by  J,  then  clearly 


and  since  (Barcilon  et  al.,1987) 


||V^l|<2^/2|j^||(l+  1  )i/2 

€  Al 


it  follows  that 


(4.1) 


R  '  e^Ai 

Thus,  as  the  Rossby  number  increases  there  is  a  tendency  for  the  vorticity  to  remain 
constant  along  streamlines. 

In  view  of  the  above  tendency,  we  look  for  a  solution  as  an  asymptotic  series  in  inverse 
powers  of  R  of  the  form 

V'  =  4'  +  +  •  •  •  1 

(4.2) 


u;  =  .^(^)+  +  ... 

K 


Also,  to  avoid  unnecessary  complications  at  this  stage,  we  restrict  ourselves  to  wind 
drivings  which  are  negative  throughout  the  entire  basin.  This  insures  the  existence  of  a 
single  gyre.  To  leading  order,  the  dynamical  equation  (1.2)  is  identically  satisfied.  To  the 
next  order,  we  get: 


9 


(4.3) 


If  ds  represents  an  element  of  arc  length  along  a  streamline,  then  after  integrating  (4.3) 
around  a  streamline,  we  deduce  that: 


Fds 


Z{<!/)=e- 


/ 


da 


(4.4) 


This  determines  the  distribution  of  vorticity  to  leading  order,  provided  that  is 
known.  Note  that  the  /3— effect  is  not  felt  to  that  order.  In  order  to  find  we  must  solve 
the  following  boundary  value  problem: 

J 


Fds 


=  €' 


IV’J'! 


in 


(4.5) 


I 

J 

^  =  0,  on  dCl 
This  boundary  value  problem  is  highly  non  standard:  not  only  is  it  nonlinear  but 

also  non  local.  Similar  problems  have  been  encountered  by  Batchelor  (1956)  and  Rhines 

&  Young  (1983)  among  others.  In  both  of  these  papers,  the  authors  were  interested  in 

determining  the  vorticity  distribution  in  a  fluid  region  with  closed  streamlines.  We  are  not 

aware  of  any  mathematical  analysis  of  such  problems.  In  section  5,  we  reformulate  this 

problem  in  terms  of  natural  coordinates  in  the  hope  that  they  will  help  future  analyses. 

For  the  time  being,  we  assume  that  $  has  been  found  and  that  it  is  unique,  and  proceed 

with  the  formal  asymptotic  expansion. 

For  convenience,  we  define  the  average  h  of  a  function  h  as  follows; 

hds 


h  = 


/ 


VVPl 


/ 


ds 

Wi 


(4.6) 


Using  now 


(4.7) 
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to  denote  the  leading  order  velocity,  we  rewrite  (4.3)  thus: 

-  {dZ/d'^)  +y}  =  F-F  (4.8) 

OS 

In  the  above  equation,  we  have  explicitly  taken  into  account  the  fact  that  F  is  negative 
and  hence  that  the  flow  around  ajiy  streamline  is  in  the  clockwise  direction  whereas  s  is 
measured  in  the  counterclockwise  one.  Also,  we  wrote  Q  for  the  magnitude  of  Q.  Clearly 

+  +C(^)  (4.9) 

Jo  Q 

In  order  to  determine  C'(^),  we  turn  to  the  second  order  equation,  viz. 

J('t,u;(2')  +  J(V>^^^u;<^^)-t- =  Q  (4.10) 


or  equivalently 


— Q  — —  {dZ / d^)ip^^^y  +  -h  y)  +  =  0 

os 


(4.11) 


Dividing  by  Q  and  integrating  around  '$'=constant  lines,  we  see  that 


^  J +  y)  ds  +  ^  Q  ds  =  0 


which  after  we  make  use  of  (4.9)  implies  that 

_  fQ-^J{¥^\S;Q-HF-F)ds) 

'  ^  fQ-Us 


jQ~H-y  +  {dz/d^W^)  +  Q-\F-  F]]  ds 

fQ-ida 


(4.12) 


We  can  now  obtain  an  equation  for  if  we  recall  that 


We  close  this  section  by  returning  to  the  two  special  forcings  discussed  previously,  namely 
the  constant  wind  curl  F  =  —1  and  the  wind  curl  equal  to  the  first  eigenfunction  of  the 
Laplacian  for  the  ocean  basin  F  =  —ui- 


ifi 


For  F  =  —1,  we  immediately  deduce  that  F  =  —1.  Therefore,  (4.5)  simplifies  greatly 


and  becomes: 


=  -e-^  in  fi, 


^  =  0,  on  dQ 


(4.13) 


For  F  =  —ui,  we  can  see  that  the  assumption  that  ^  is  proportioned  to  does  not 


lead  to  any  inconssitencies.  Indeed,  under  this  assumption,  the  problem  for  'i'  reduces  to 


=  -e  in  Q, 


^  =  0,  on  dCl 


(4.14) 


and  consequently 


(4.15) 


5.  The  problem  for  We  return  to  the  problem  for  the  leading  term  in  the 


eisymptotic  series,  namely 


=  F  in  n 


=  0  on 


We  restrict  our  attention  to  forcings  of  one  sign,  say  F  <  0.  To  simplify  the  right  hand  side. 


we  switch  to  the  natural  coordinates  ^  and  s,  where  s  is  the  arc  length.  The  unknowns  are 


X(3,  ^)  and  T(3,  ^)  which  provide  a  parametric  representation  of  the  streamlines.  Define 


the  Jacobian 


d{X,Y) 


-.V.'  '.■'•.■V 


?5!l 


SfM 


' ) 


1 


m 


cl 


#  ii 


For  future  reference,  we  record  the  following  formulas: 


=  -r'^Y, 


=  /-'A', 


1 5,  =  r'^Y^ 

1  5,  = 


As  a  result,  we  can  deduce  the  expression  for  the  Laplacian  of  ^ 


■1  (  Y  — 


or  better  still 


^  -  (A,A*  ^  -  (A»A»  + 


In  the  process  of  deriving  the  above  expression,  we  have  made  use  of  the  fact  that 


A^  +  y/  =  1 


From  the  definition  (4.1)  of  J,  we  have 


I  ^  X,Y^  -  X^Y, 


Finally,  we  note  that 


Q  =  jn  +  n  =  i' 


As  a  result,  the  problem  (5.1)  is  transformed  into  the  solution  of  the  three  coupled  partial 


differential  equations 


^  -(A,A^  +  y,F^)^  -(A>A..  +  y.,y„)^  =  e-'i^^, 


I  =  X,Y<s.  -  X<iY,, 


A^  +  y,^  =  1. 


6.  Discussion:  We  have  seen  that  the  lack  of  a  priori  bound  on  the  L°°  norm  of 


S>!l 


the  gradient  of  vorticity  prevented  us  from  obtcuning  a  uniqueness  result  veJid  beyond  the 
viscous,  small  R  regime.  At  this  stage,  we  do  not  know  whether  such  a  bound  is  just 
difficult  to  derive  or  simple  impossible,  i.e.  whether  one  can  envisage  wind  curl  which  lead 
to  very  large  values  of  the  gradient  of  vorticity. 

The  same  remarks  apply  to  the  /—plane  analog.  However,  in  this  case,  we  saw  that 
for  flows  such  a.s  those  driven  by  a  wind  curl  which  is  either  constant  or  proportional  to 
the  first  eigenfunction  of  the  basin,  the  vorticity  gradient  is  bounded  and  the  solution 
is  unique  for  all  values  of  the  Rossby  number.  Incidentally,  for  these  two  drivings  the 
majdmum  value  of  the  vorticity  was  in  the  ocean  rather  than  along  the  coast  as  was  the 
case  in  the  numerical  calculations  of  Veronis  (1966b)  for  the  /3— plane. 


For  large  values  of  R,  we  have  derived  a  non-local  problem  for  the  leading  order 
strearafunction.  The  vorticity  is  then  roughly  constant  on  these  streamlines.  The  associ¬ 
ated  circulation  shows  a  very  weak  East- West  eisymmetry. 


>**?s*!i 
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